Proceedings of the 6th European VLBI Network Symposium 
Ros, E., Porcas R.W., & Zensus, J. A. (eds.) 
June 25th-28th 2002, Bonn, Germany 




EUROPEAN 
NETWORK 



General relativistic model for experimental measurement of the 
speed of propagation of gravity by VLBI 

S. Kopeikin 1,3 and E. Fomalont 2 



O 
O 

(N 

C 



> 

(N 
(N 
O 

o 
o 

cr : 



1 Department of Physics and Astronomy, University of Missouri - Columbia, 223 Physics Bldg., Columbia, Missouri, 
65211, USA 

2 National Radio Astronomical Observatory, 520 Edgemont Road, Charlottesville, VA 22903, USA 

3 E-mail: kopeikins@missouri.edu 

Abstract. A relativistic sub-picosecond model of gravitational time delay in radio astronomical observations is 
worked out and a new experimental test of general relativity is discussed in which the effect of retardation of 
gravity associated with its finite speed can be observed. As a consequence, the speed of gravity can be measured 
by differential VLBI observations. Retardation in propagation of gravity is a central part of the Einstein theory 
of general relativity which has not been tested directly so far. The idea of the proposed gravitational experiment 
is based on the fact that gravity in general relativity propagates with finite speed so that the deflection of light 
caused by the body must be sensitive to the ratio of the body's velocity to the speed of gravity. The interferometric 
experiment can be performed, for example, during the very close angular passage of a quasar by Jupiter. Due to 
the finite speed of gravity and orbital motion of Jupiter, the variation in its gravitational field reaches observer on 
Earth not instantaneously but at the retarded instant of time and should appear as a velocity-dependent excess time 
delay in addition to the well-known Shapiro delay, caused by the static part of the Jupiter's gravitational field. Such 
Jupiter-QSO encounter events take place once in a decade. The next such event will occur on September 8, 2002 
when Jupiter will pass by quasar J0842+1835 at the angular distance 3.7'. If radio interferometric measurement of 
the quasar coordinates in the sky are done with the precision of a few picoseconds (~ 5 /xas) the effect of retardation 
of gravity and its speed of propagation may be measured with an accuracy about 10%. 



1. Theoretical Background 

Experimental verifications of the basic principles un- 
derlying Einstein's general relativity theory are impor- 
tant for fundamental physics. All previous experimen- 
tal tests of general relativity in the solar system have 
relied upon the static Schwarzchild solution (Will 1993) 
and, therefore, were not sensitive to the effects entirely 
associated with the propagation speed of gravity. It is 
worth noting that gravitational waves are inherent to 
the radiative (far) zone of a system emitting the waves 
flMisner, Thornc fc Wheeler 197$ parish fc Weiss 1999| ). 



However, the gravitational waves do not propagate freely 
through the interior of a non-radiative (near) zone of the 
system. Nevertheless, the process of generation of gravi- 
tational waves produces retarded effects in the near zone 
leading to appearance of the gravitational radiation re- 
action force in the relativistic equations of motion of ex- 
tended bodies comprising a self-gravitating astronomical 



veal that electromagnetic signals interact with the light- 
ray deflecting bodies only through the retarded gravita- 
tional fields - the observational effect which must be ac- 
counted for in precise data processing algorithms adopted 
for the microarcsecond space astrometry. 

This gravitational light-propagation theory, in the case 
of a static spherically-symmetric field, gives the same re- 
sult as that predicted by Einstein for the bending of light, 
as — 4:GM/(c 2 R6), where M is the mass of the light- 
ray deflecting body, R is the distance from observer to 
the body, and 9 is the (small) angle in the sky between 
the undisturbed geometric positions of the source of light 
and the center of mass of the massive body. Furthermore, 
this theory allows to calculate the correction, ape, to 
the Einstein deflection ole related to the variability of 
the gravitational field produced by the motion of the 
light-ray deflecting body. We were successful in proving 



system (Damour et al. 1989). Existence of this force is a 



consequence of the finite speed of propagation of gravity 
as it was experimentally confirmed by Taylor (1994). 



We have found (Kopeikin 2001) that the gravitational 



bending of light passing through the gravitational field of 
a moving massive object like Jupiter, though being domi- 
nated by the spherically-symmetric component of its grav- 
itational field, also contains terms associated with the fi- 
nite speed of propagation of gravity. Our calculations re- 



( Kopeikin 2001) that these corrections in the bending of 
light are inherently associated with the finite speed of 
propagation of gravity and in case of slowly moving bod- 
ies can be parameterized as apo ~ (1 + 5){cie / 0){v / c), 
where v is the orbital velocity of the light-ray deflecting 
body with respect to the barycenter of the solar system 
projected on the plane of the sky, and <5 = c g /c — 1 is a 
fitting parameter used in data analysis. It is chosen such 
that 6 = 0, if the speed of gravity c g equals the speed of 
light c. 
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Parameter S is a close analogue of the parameter 
a 2 = ( c g/ c ) 2 — 1 °f the parameterized post-newtonian 
(PPN) formalism which quantifies possible violation of 

the local Lorentz invariance ( |Will 1993Q F1- K 

was shown 



([Nordtvedt 1987|) that a 2 <4x 10^ under the (rather 



restrictive) assumption that the preferred frame is real- 
ized by the cosmological Hubble flow. If one abandons any 
anthropic assumption about the speed of the solar system 
with respect to an (actually unknown) preferred frame the 
limit ona2 < 0.1 can be obtained from the an alysis of th e 
anomalous perihelion shifts of inner planets ([Will 1993j) . 
The primary purpose of our experiment is, however, to 
observe directly the effect of retardation in propagation of 
gravitational field in the solar system rather than improv- 
ing limits on ai- Nevertheless, we would like to empha- 
size that relativistic effects in propagation of light through 
time dependent gravitational fileds are also sensitive to vi- 
olation of local Lorentz invariance. 

The largest measurable contribution to the variable, 
time-dependent part of the solar system gravitational field 
comes out from the orbital motion of Jupiter. The min- 
imal value of the impact parameter of an incoming light 
ray from a quasar that can be achieved for Jupiter is also 
the least possible amongst all the solar system bodies. 
Therefore, it is highly sensible to make an attempt for 
detection of the "gravity retardation" effect by observing 
very accurately the gravitational deflection of light from 
a background source (quasar) caused by the motion of 
Jupiter aro und the baryce nter of the solar system. As ex- 
plained in ( Kopeikin 2001 ) the magnitude of the observed 
effect is directly translated to the measured value of the 
propagation speed of gravity c g . This is the essence of the 
new test of general relativity which has never been done 
before with sufficient accuracy. 

Radio astronomical methods of VLBI are the most 
accurate for measuring gravitational deflection of elec- 
tromagnetic waves. The two most precise measure- 
ments of the bending of radio waves near the sun 



(Lebach et al. 1995; Robertson, Carter & Dillingcr 1991) 
were accurate to about 0.1%. However, these observations 
were insensitive to the speed of gravity effect ape because 
of the relatively large impact parameter of the incoming 
light ray. Even at the solar limb, the magnitude of ape 
is ~ 10~ 5 of the static gravitational bending = 1.75' 
and is totally unobservable because of the highly turbulent 
solar magnetosphere. 

On September 8, 2002 Jupiter will pass at an angu- 
lar distance of 3.7' from the quasar J0842+1835 making 
an ideal celestial configuration for measuring the speed 
of propagation of gravity by using the phase-referencing 
VLBI technique ( [Fomalont fc Kopeikin 2002j ). The en- 
counter in 2002 is especially favorable because: (1) it oc- 
curs when Jupiter is relatively far from the Sun (the next 
near occultation which occurs is a few degrees from the 
Sun), and (2) the five critical hours of the closest approach 



occur when Jupiter is near the transit line for VLBA ob- 
servations. 

We have estimated that for this Jupiter-quasar en- 
counter the deflection from the static gravitational field 
and from the propagation of gravity are, respectively, 
ctE — 1-26 mas and apg = 53 /ias ^ both in the plane 
of the sky with the static bending radially from Jupiter 
and the propagation bending in the direction of Jupiter's 
motion (see Eqs. (0) and (|J) in Sec. 3). As one can see the 
ratio \apc\/oLE — 0.04 is much larger for Jupiter than for 
the Sun which is explained by the ability to get a smaller 
impact parameter 9 for the light ray passing by Jupiter 
than that for the Sun. 



2. Relativistic Model of VLBI Time Delay 

Detection of the effect of gravity propagation requires 
a more advanced VLBI model for light propagating in 
the time dependent gravitational field of the solar sys- 
tem. Such a model must be valid to a precision of bet- 
ter than 0.1 ps. In the present paper we discuss the ap- 
propriate corrections to the standard Shapiro time delay 
(Shapiro 1967) which bring the accuracy of the model up 
to the necessary threshold. 

The general formula for the relativistic time delay AT 
in the field of a system of moving bodies is given in 
(|Kopeikin 200l|) 



AT 



G N 

a=l 



;i-ifc-M0) T(Qrfc 



, (1) 



where T(£) = 1/^/1 — c 2 t"a(C) is the Lorentz factor, 



is the PPN parameter ( Will 1993| ), m a is the mass of the 
ath body, t* is the time of the closest approach of electro- 
magnetic signal to the barycenter of the Solar system PI 
x a (t) are coordinates of the ath body, v a (t) = dx a {t)/dt is 
the (non-constant) velocity of the ath light-ray deflecting 
body, k is the unit vector from the point of emission to 
the point of observation, s is a retarded time obtained by 
solving the gravitational null cone equation for the time 
of observation of photon t = s + cZ 1 \x — x a (s)\, and Sq is 
found by solving the same equation written down for the 
time of emission of the photon to = sq + c~ 1 \xq — x a (so)\. 

One emphasizes that the equations for the retarded 
times depend on the speed of gravity c g , but not the speed 
of light c. This is because they were obtained by solving 
Einstein equations for the space-time metric perturbations 
by m aking use of retarded Liena rd-Wiechert tensor poten- 
tials ( Kopeikin fc Schaier 1999 ). These retarded gravita- 
tional potentials describe propagation of gravity without 
any relation to the problem of propagation of light in the 
gravitational field. In general relativity c g = c numerically. 



One notices that 8 — 012/2 in the first approximation. 



2 It corresponds to the time delays 122.2 and 5.1 picoseconds 
respectively on a baseline b = 6000 km. 

3 The time t* is used in calculations as a mathematical tool 
only. It has no real physical meaning because of its dependence 
on the choice of a coordinate system. 
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However, when light propagates through time-dependent 
gravitation field, in principle, one can separate relativistic 
effects associated with propagation of light and gravity. 

The Earth and Sun also contribute significantly to the 
gravitational time delay and must be included in the data 
processing algorithm in order to extract accurately the 
effect of retardation of gravity. Precise calculation of the 
integral (JTj) for two radio antennas gives the differential 
VLBI time delay 



A 2 T- AiT 



A + A, 7 + Ajp G 



(2) 



The first term in the right hand side of ^ describes the 
gravitational (Shapiro) time delay due to the gravitational 
field of the Earth 

GM e Xi + K ■ X\ 
A e = (I+7) — — In . 

c J X 2 + K ■ A 2 



It can reach 21 ps for the baseline b — 6000 km. 

The second term in the right hand side of (|J) describes 
the gravitational (Shapiro) time delay due to the Sun 



A© = (1+7) — 5— In 



K r 



1© 



r 2Q + K ■ r 2G 



(4) 



It can vary (for b = 6000 km) from 17 x 10 4 ps for the light 
ray grazing the Sun's limb to only 17 ps when direction 
to the source of light is opposite to the Sun. 

The third term in the right hand side of (||) is the 
Shapiro time delay due to the static part of the gravita- 
tional field of Jupiter 



A., 



GMj , r u + K -r u 
(I+7) — — (1 + K ■ vj) In z ; — . (5) 



r 2 j + K r 



2.1 



Finally, the forth term in the right hand side of (||) 
is the time delay caused by the finite speed of gravity as 
predicted by general relativity theory ( Kopeikin 2001 ) 



A JPG = 2(1 + 6) — — . 6 

c ru + K ■ ru 

In formulas we use the following notations: M® 

- mass of the Earth, M Q - mass of the Sun, Mj - mass 
of Jupiter, Vj(ti) - the barycentric velocity of Jupiter, 
and K - the unit vector from the barycenter of the so- 
lar system to the quasar observed. Also, for each i = 1, 2, 
one has the baseline vector b = Xi — X 2 , r iQ = |r iQ |, 
n.j = \ru\, Nu = ru/ru, r iG = Xi(t{) - x G (U), 
r u = Xi(ti) - xj(U), x. t = x e (tj) + Xi(U), where 
Xi(ti) are the geocentric coordinates of i-th VLBI sta- 
tion, x^ - the barycentric coordinates of the geocenter, 
Xq - barycentric coordinates of the Sun, xj - barycen- 
tric coordinates of Jupiter, and ti is time of arrival of the 
plane front of electromagnetic wave from quasar to the ith 
VLBI station. 

We notice there are two relativistic parameters to be 
measured in order to test validity of general relativity the- 
ory - the PPN parameter 7, and the speed of gravity pa- 
rameter S = Cg/c— 1. The best experimental measurement 
of parameter 7 had been conducted by Lebach et al. (1995) 



who obtained 7 = 0.9996 ± 0.0017 in an excellent agree- 
ment with general relativity. The primary goal of the new 
experimental test of general relativity is to measure the 
parameter S which will set up limits on the numerical value 



of the speed of gravity c g (Fomalont & Kopeikin 2002). 

During the passage of Jupiter near the quasar the time- 
dependent impact parameter £(£) of the light ray with re- 
spect to Jupiter will be always small as compared with 
the distance from Earth to Jupiter which will be approx- 
imately 6 AU. It is convenient to introduce the unit vec- 
tor n = £/|£| along the direction of the impact parame- 
ter according to definition sin 9 n = (K x (Nu x K)), 
where 9 is a small angle between the unperturbed as- 
trometric position of the quasar and that of Jupiter. 
Making use of the (impact parameter 9) expansion Nu = 
(3) -(l-^y^K + fln+O^ 3 ), we obtain the functional struc- 
ture of the Shapiro time delay A j and the speed of gravity 
delay Ajpq in a more explicit form (assuming for simplic- 
ity 7=1) 



A., 



AGMj 
c 3 ru 



n B {n-Bf (K x B) 2 



e 2 



2ru6 



(7) 



A JPG = (l + S)^ b -^-^ K - b \ (8) 



c 4 ru 



where B = b c~ l (K ■ b)(v 2 - v,j) + 0{c~ 2 ) , and all 
quantities in the right sides of Eqs. (0)-(||) are taken at 
the time t\. 



3. The Effect of the Magnetosphere of Jupiter 

In addition to various special and general relativistic ef- 
fects in the time of propagation of electromagnetic waves 
from the quasar to the VLBI antenna network, we must 
account for the effects produced by the jovian magneto- 
sphere. Measurements obtained during the occultations of 



Galileo by Jupiter indicate (Flasar et al. 1997) that near 
the surface of Jupiter the electron plasma density reaches 
the peak intensity No = 1.0 x 10 10 m~ 3 . We shall assume 
that the jovian magnetosphere is spherical ^] and a ra- 
dial drop-off of the plasma density N(r) is proportional to 
l/r 2+A where r is the distance from the center of Jupiter. 
The guess is that A > 0, and we will assume that A = for 
the worst possible case. Hence, radial dependence of the 
electron plasma density is taken as N(r) — Nq(Rj /r) 2+A , 
where Rj = 7.1 x 10 7 m is the mean radius of Jupiter. 

The plasma produces a delay AT in the time of prop- 
agation of radio signal which is proportional to the col- 
umn plasma density in the line of sight given by integral 
( |Yakovlev 198^ ) 



Ni = 



N(r) dr 



n A-\- \ ^ ^,2 



d 2 



N(r) dr 



(9) 



4 In reality the magnetosphere has a dipole structure and 
we speculate that our model which assumes circular symme- 
try grossly underestimates the plasma content along the polar 
direction where the closest approach occurs. 
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where r$ and r\ are radial distances of quasar and radio 
antenna from Jupiter respectively, and d = |£| is the im- 
pact parameter of the light ray from the quasar to Jupiter 
0. In the experiment under discussion the impact param- 
eter is much less than both rp and r%. Hence, 



^(m- 2 ) = N Rj 



Rj_ 
d 



r(i + 4) 



(10) 



where T(z) is the Euler gamma-function. The plasma time 
delay 



AT(s) = 40.4c~V~ 2 iV ; 



(11) 



where c is the speed of light in vacuum measured in m/sec, 
v is the frequency of electromagnetic signal measured in 
Hz. 

It is worthwhile noting that, in fact, the VLBI array 
measures difference in path length between the radio tele- 
scopes. Hence, one has to differentiate Ni in expression 
( |TT| ) with respect to the impact parameter d and project 
the result on the plane of the sky. This gives a magneto- 
spheric VLBI time delay of 

A, M (ps) =6.3xl0-^ + l)§(^) 2 ^, (12) 

d \ v / c 

normalized to the frequency vq — 8.0 GHz. Substituting 
d = 13Rj and taking the baseline b — 6000 km we find 



A JM = 2.34 (V^f ps (A = 0) 

A JM = 0.13 (V) 2 PS (A = 1) 

A JM = 0.03 [u /uf ps (A = 2) 



(13) 
(14) 
(15) 



This represents the pure bending from Jupiter's magne- 
tosphere which should be compared with the propagation 
of gravity time delay Ajpg = 5.1 ps at closest approach. 

If we observe at the dual frequencies at 2.3 GHz and 
8.4 GHz in the normal geodetic mode, we certainly can 
determine the ionospheric (both from Jupiter and from 
the Earth) effects. However, our sensitivity at 8.4 GHz 
will be decreased because of only one polarization and 
only half the total bandwidth. The noise in the iono- 
sphere/magnetosphere bending may also be a limit. As a 
rough order of magnitude, the position error per day is 10 
/xas at 8.4 GHz and 30 /las at 2.3 GHz. Whatever bending 
we obtain at 2.3 GHz, about 10% is removed from the 8.4 
GHz bending. Thus, the ionosphere/magnetosphere cor- 
rection will have an error of 3 to 4 ^as. 

Only in the worst case scenario will the effects of 
Jupiter's magnetosphere be significant at 8 GHz observ- 
ing. Perhaps we should observe at two widely spaced fre- 
quencies in the 8 GHz band, say at 8.0 GHz and 8.5 GHz, 
then, independently reduce the data at the two frequen- 
cies. The gravitational bending delay (]?]) and gravitational 
retardation of gravity delay (0) are both independent of 
frequency. Any plasma delay should scale inversely with 

5 This impact parameter d ~ 14_Rj on September 8, 2002 
and corresponds to the angle 6 — 3.7' between the quasar and 
Jupiter in the plane of the sky. 



the frequency-squared and can be determined by looking 
at the difference measurements. However, this small fre- 
quency difference will produce very large errors in the es- 
timate of the jovian bending and also be strongly affected 
by Earth ionospheric contamination. 

If we observe at 15 GHz instead, the magnetospheric 
and ionospheric delays are both a factor of four smaller 
and almost certainly negligible. However, the system sen- 
sitivity is less and one of the calibrators may be too weak 
to reliable detect (J0839+1802). We are still unsure about 
the most optimum method to deal with the possible jovian 
magnetosphere component. 

Time variability of the Jupiter magnetosphere could 
cause problems. For example, if there were very large, 
chaotic changes in the jovian magnetosphere, then we 
could lose coherence over a minute of time. However, this 
fluctuation model is very pessimistic and unlikely, and 
would probably average out to the steady state model. 

4. Summary 

We believe that the differential VLBI experiment in 
September 2002 can measure the retardation effect in 
propagation of gravity and determine the speed c g of its 
propagation with 10% to 20% accuracy. If the experiment 
is successful it will provide a new independent test of gen- 
eral relativity in the solar system. 
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